PHASE TRANSITIONS FOR THE BRUSSELATOR MODEL 



TIAN MA AND SHOUHONG WANG 

Abstract. Dynamic phase transitions of the Brusselator model is carefully 
analyzed, leading to a rigorous characterization of the types and structure of 
the phase transitions of the model from basic homogeneous states. The study 
is based on the dynamic transition theory developed recently by the authors. 



1. Introduction 

Bclousov-Zhabotinsky (BZ) reactions are now one of a class of reactions that 
serve as a classical example of non-equilibrium thermodynamics, resulting in the 
establishment of a nonlinear chemical oscillator. 

The main objective of this article is to study the dynamic phase transitions of 
the Belousov-Zhabotinsky reactions focusing on the Brusselator, first introduced by 
The Brusselator is one of the simplest models in nonlinear chemical systems. 
It has six components, four of which retain constants, and the other two permit 
their concentrations vary with time and space. The chemical reaction consists of 
four irreversible steps, given by 

B + X % Y + D, 
2X + Y 3X, 
X^>E. 

where A and B are constant components, D and E are products, and X and Y are 
the two components variable in time and space. Over the years, there have been 
extensive studies for the Brusselator and related chemical reaction problems; see 
among many others [T3J [21 HI SI H2] and the references therein. 

In this article, we address the dynamic phase transition of the Brusselator model. 
In particular, we derive a complete characterization of the transition from the ho- 
mogeneous state. There are two aspects of this characterization. First our analysis 
shows that both the transitions to multiple equilibria and to time-periodic solutions 
(spatiotemporal oscillations) can occur for the Brusselator model, and are precisely 
determined by the sign of an explicit nondimensional parameter So — 5% as defined 
by (33) and (jXtJl) . 

Then in both transition cases, the dynamic behavior of the transition is classi- 
fied based on the new dynamical classification scheme, introduced as part of the 
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dynamical transition theory developed recently by the authors; see [SJ|6l[7]. With 
this classification scheme, phase transitions are classified into three types: Type-I, 
Type-II and Type-Ill, which, in more mathematically intuitive terms, are called 
continuous, jump and mixed transitions respectively. For the Brusselator, the dis- 
tinction of the dynamic transition types are determined, again, by the signs of some 
nondimensional computable parameters. 

It is worth mentioning that the main philosophy of the dynamic transition theory 
is to search for the full set of transition states, giving a complete characterization 
on stability and transition. The set of transition states is represented by a local 
attractor. Following this philosophy, the dynamic transition theory is developed 
to identify the transition states and to classify them both dynamically and phys- 
ically. With this theory, many long standing phase transition problems are cither 
solved or become more accessible, providing new insights to both theoretical and 
experimental studies for the underlying physical problems. 

This article is organized as follows. Section 2 introduces the model and its 
mathematical set-up, and Section 3 addresses the principle of exchange of stabili- 
ties. Dynamic transitions of the model are addressed in Sections 4-6, with physical 
remarks of the main results given in Section 7. 



Let u\,U2,a and b stand for the concentrations of X,Y,A and B. Then the 
reaction equations of (|l.ip read 



2. The Model and its Mathematical Set-up 



(2.1) 




To get the nondimensional form of (|2.ip . let 



t — t 7 

a = (kl/klk 3 y/ 2 a, 




b = (k 4 /k 2 )X, 



a l = fk^. 



for i = 1,2. Omitting the primes, the equations (|2.1[) become 




^1,^2, a, A > 0. 



The equations ()2.2[) have a constant steady state solution 

uq = (a, A/a). 



Make the translation 



u 1 =a + Vt, 



A 

u% = - + v 2 , 



a 
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then the equations (|2.2|) are written as 

dv\ , , 2 2A 2 2 
= /xi A«i + (A — l)vi + a v-i H + 2ctoiv 2 + ^^2, 

(2.3) 



dt n 



9V2 A \ 2 2A 2 2 

at a 

There are two types of physically-sound boundary conditions: the Dirichlet 
boundary condition 

(2.4) v= (vi,v 2 )=0 on<90, 
and the Neumann boundary condition 

(2.5) — = on dQ. 

On 

Define the function spaces 

H = L 2 (n,R 2 ), 

( H 2 (n,R 2 )DH^(n,R 2 ) for b.c flU), 

Hi = < g w 

[{«G JJ 2 (fi,R 2 )| — |aa = 0} for b.c. (1231). 

Define the operators La = A + B\ and G : H± H by 

= (/jiAt>i,/i 2 Au2), 
^ 2 ^ -B A w = ((A - + a 2 v 2 , -\vi - a 2 v 2 ) 1 

— + 2auiW2 + «i — 2awiW2 — V1V2 

Thus the equations (|2.3j) with (|2.4p or with (|2.5j) can be written in the following 
abstract form 

dv 

(2.7) — =L x v + G(v,X). 

3. Principle of Exchange of Stability (PES) 
Consider the eigenvalue quations of (|2.3p 
^ ^ fJ-iAvi + (A - + a 2 u 2 = fax, 

/j, 2 Av2 — Aui — a 2 v 2 — fiv2, 



with the boundary condition (I2.4[) or (I2.5[) . 

Let pfc and efc be the kth eigenvalue and eigenvector of the Laplacian with either 
the Dirichlet or the Neumann condition: 

Ae k = -p k e k , 

(3-2) de k 

e k \an = or — — \ m = 0. 

On 

Denote by M k the matrix given by 



-/Ltipfc + A - 1 

-A - ct z 



M k = [ .. : .2 , fc=i,2,... 
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It is clear that all eigenvalues and eigenvectors (j>j: of (|3.1[) satisfy the following 
equations 

(3.3) ^ = / efc ; ± 

where G K 2 are the eigenvectors of M k , (3^ are the eigenvalues of M k , which are 
expressed as 

(3-4) (A) =i[A - (/iiPfe + M2Pfe + « 2 + 1)] 



It is clear that j3 k (A) < /3^(A) = if and only if 



± - [(A - mpk - p-2pk - ct 2 - l) 2 
+ <i(Xp 2 pk - (pipk + l)(M2Pfe + a 2 ))] 



1/2 



A = — (/iipfc + l){p%Ph + a 2 ), 

A < ^iip fe + /i 2 pfe + a 2 + 1, 
and /3^(A) = ±crfe(A)i with er^. ^ if and only if 

A = HiPk + P-2Pk + a 2 + 1, 

A < (MiPfc + l)(M2Pfe + a 2 )- 

Thus we introduce two critical numbers 

(3.5) A = min {p\Pk + ^){P2Pk + a 2 ), 

Pk P2Pk 

(3.6) Ai = pipi + p 2 pi + a 2 + 1. 
Obviously, the following lemma holds true. 

Lemma 3.1. Let Xq and Ai be the two numbers given by \3. 5\) and \3. 6\) . Then we 
have the following assertions: 

(1) Let Xq < Ai, and ko > 1 be the integer such that the minimum is achieved 
at pk in the definition of Xq . Then /3^~ (A) is the first real eigenvalue of 
\3.1\) near X = Xq satisfying that 

( < if X < X 
f3+(X)\ = ifX = X VfcGN with p k = p ko , 
( 3J ) [ > ifX > A 

Repf (A ) < Wfjf f (3+ with Pk = p ko . 

(2) Let Ai < Ao- Then /3^~ (A) = /3-f (A) are a pair of first complex eigenvalues 
of \3. 1\) near X = X\, and 

( <0 i/A<Ai, 
Re/3+(\) = Refc(X){ = ifX = X u 
( 3 - 8 ) [ >0 ifX>X u 

Re^(X 1 )<0 Vfc>l. 

Remark 3.1. (A) are simple complex eigenvalues at Ai(< Ao), and in general, 
if p ko is a simple eigenvalue of (|3.2[) . then (3 ko (X) is also simple at Ao(< Ai). 
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4. Transition from real eigenvalues 

Hereafter, we always assume that the eigenvalue /3~£ in (|3.7j) is simple. Based 
on Lemma 13.11 as Ao < Ai the transition of (I2.7|) occurs at A = Ao, which is from 
real eigenvalues. Let pk be as in Lemma 13 - 1 1 and ek the eigenvector of 
corresponding to efc satisfying 



(4.1) / ei o dx^0. 

Then, under the condition (I4.ip . for the system (|2.3[) with (|2.4I) or with (I2.5[) we 
have the following transition theorem. 

Theorem 4.1. Let Xq < X±. Then the system \2. 7[ ) has a transition at X = A , 
which is mixed (Type-Ill). In particular, the system bifurcates on each side ofX = A 
to a unique branch v x of steady state solutions, such that the following assertions 
hold true: 

(1) On X < Xq, the bifurcated solution v x is a saddle, and the stable manifold 
T\ of v x separates the space H into two disjoint open sets and U^, such 
that 9 = 06 Ui is an attractor, and the orbits of inU^ are far from 
v = 0. 

(2) On X > Xq, the stable manifold T° of v = separates the neighborhood O 
of u = into two disjoint open sets 0\ and 0\, such that the transition is 
jump in 0\, and is continuous in 0\- The bifurcated solution v x € 0\ is 
an attractor such that for any ip € 0\ , 



lim \\v(t,tp)-v*\\ H = 0, 



where v(t, ip) is the solution of with v(0, ip) = tp. 
(3) The bifurcated solution v x can be expressed as 



v X = CPl(X)Z+e k0 +o(\f3t \), 
(-M2/?fc .MlPfeo + !)j 

(o^2PfcoG»2Pfc + a 2 ) - Q^iPfcp + 1)) fn e 2 ko dx 



/ 42 j Cko = (-M2pfco.MlPfeo + 1), 



c 



Proof. We apply Theorem A. 2 in [S] to prove this theorem. Let $ be the center 
manifold function of (I2.7[) at A = A . We need to simplify the following expression: 

(4-3) g(3,) = 1 (G(yri o +Hy)),&i 

where y e M 1 , G is the operator defined by (|2.6j) . (j)^ g is the eigenvector of ()2.5[) 
corresponding to P^ (Xq) = 0, and (j>^* is the conjugate eigenvector. By (|3.3|) . 

i) « 2 V £ 



(4.4) 


n = £,t e k , 


(4.5) 


( Aq - (^ip fel 

V -Ao 


(4.6) 


/ Ao - (a*iP/c, 
I a 2 



fco i 



-(M2Pfc +" 2 ) / V C 2 



= 0. 
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By definition of Ao and fco , we infer from (|4.5|) and (|4. 6[) that 

(4-7) tl = (# ol ,# o2 ) = (-»2Pk ,»iPk + 1), 

(4-8) = (Hi\,ti; 2 ) = (M2P.„ + a 2 , a 2 ). 

Denote by o{k) — o(\y\ k ). By $(y) = o(l), the function g(y) in (|4.3[) is rewritten as 

9(V)= (6+ l ^ (G{yrt o Ui;)+o{2). 
By (gU) and dUTJ) we see that 

2y2 ("o^ 2 ^ ~ a ^P*o{l^Pk + 1)J e+fco + o(2), 



2y2 ^ ^^^fe _ a(i 2 Pk (viPko + 1) ) ejfc 4 ,H2). 



Thus we deduce from (l4~4l) and (t4~7ll - (l4~8l) that 



(^fco'^fco*) = ( a2 (ViPk„ + 1) - P2Pk {P2Pk + a 2 ) / e^cfe 
(G(#+, 0+*) = y 2 ■ ^L( MlPfco + 1) / e | o dr + o(2). 



Q 

Therefore the function f|4.3[) is given by as 



Q 



.9(2/) = -V + °(2), 



c 



and the theorem follows from Theorem A. 2 in [8]. The proof is complete. □ 

Remark 4.1. If the domain Vt ^ (O.L)xD with DcR 2 being a bounded open set, 
then the condition (|4.ip holds true for both the Dirichlet and Neumann boundary 
conditions (|2.4[) and (12.51) . If Q = (0, L) x D, then (|4.ip is not true for the Neumann 
condition (|2.5[) . and is not true for the Dirichlet condition (|2.4p if the number m in 
Pk = "Y?" + p' ko is even. Here p^, is an eigenvalue of the equation on D: 



— Ae = p' kg e for x E D, 
e\dD = 0. 

□ 

Now, we consider the case where ()4.1j) is not true, i.e., 

(4.9) f e 3 ko dx = 0. 

As mentioned in Remark 14.11 the condition (|4.9[) may hold true if f2 = (0,L) x D. 
We introduce the following parameter 

(4.10) b x =[oiplpl Q {niPk + 1) / el dx -2a 2 p J2 pk ^2e\ dx 

Jn Jn 

- 2(A«iPfc + l)(2M2Pfc + " 2 ) / ^iefc rfx] 



x [a 2 (^i ( o fco + 1) - M2Pfco(M2Pfc n + a 2 )] \ 
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where ijj — (ipi,ip2) satisfies 

W AV»i + (A - l)^i + o?^h = -^^(jiipko + l)e 2 ko , 

(4-H) M2 A^ 2 - Ao^ - a 2 i> 2 = 2 J^L{^ lPk0 + l)^ 

a 

dip 

Won = (or — jan = 0). 

By the Fredholm Alternative Theorem, under the condition (I4.9p . the equation 
(|4.11jl has a unique solution. 

Theorem 4.2. Let hold true, Ao < Ai, and b\ is the number given by J^. 10ty . 
Then the transition of \2. 7\ ) at X = Xq is continuous ifb\ < 0, and is jump ifb\ > 0. 
Moreover, the following assertions hold true: 

(1) If b\ > 0, \2. 7\ ) has no bifurcation on X > Xq, and has exact two bifurcated 
solutions and v_ which are saddles. Moreover, the stable manifolds 
and of the two bifurcated solutions divide H into three disjoint open sets 
U$, U$, Ut such that v = e is an attractor, and the orbits in U± are 
far from v = 0. 

(2) If b\ < 0, has no bifurcation on X < Xq, and has exact two bifurcated 
solutions vi_ and , which are attractors. In addition, there is a neighbor- 
hood O C H of v = 0, such that the stable manifold V of v — divides O 
into two disjoint open sets 0+ and such that v± £ 0+, v_ & 0_, and 
v§. attracts <D$_; 

(3) The bifurcated solutions v± can be expressed as 

4 = ±C( / 8+(A)) 1 /^+ efeo + o (| /0 +|V2 )) 

(4.12) 



c = 



-1 1/2 

where b\ is as in J^._?0[ ). 

Proof. We use Theorem A.l in [8] to prove this theorem. To get the function g(y) 
in (14.3[) . we need to calculate the center manifold function Q(y). By (A. 10) in [9], 
i>(y) satisfies 

(4-13) L Ao $ = -P 2 G(#+), 

where P 2 : H — > £2 is the canonical projection, L\ is as in (|2.6j) . 0jT o and are 
given by (|4.4j) . and 

£ 2 = {«e#| («,O = >- 

We see that 

+ f 2a"Va/4 ^ + + °( 2 )- 

ly</,fcJ ~ I - 2a-V^ (MiP fe0 + l)2/ 2 4 + 0(2). 

Let 

(4.14) $ = VV + o(2), V = (V>1,$0 G 
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By (HU), (e 2 ko ,-e 2 ko ) e E 2 . Hence, it follows from (|4J3]l and (j4H| that 

(4-15) L x Jj = -2a- 1 fi 2 2 plMPk + ^(efco.-O. 

which is an equivalent form of (|4.11l) . 
By ()4.14p . we have 

G(#+ + *) = GQ/0+ + y 2 ^) + o(3), 

' °{y 2 ) e l + y 3 [P-lpl )(PiPk + l)e| - 2ap 2 p ko ^ 2 e ko 

- 2a~ 1 Tp 1 e ko (p 1 p ko + l)(2/i 2 pfe D + a 2 )] + o(3), 
O(y 2 ) e t - y 3 [plPk (PiPko + 1)4 ~ ^ap 2 p ko tp 2 e ko 

- 2a~ 1 ipie ko (pipk + l)(2^2Pfc + a 2 )] + o(3). 



G(#+ +yV)= ^ 



Hence we deduce from fO}, P~g|) and (|43| that 



(G(^+ +$),$+*) =^^- a/^^^ + l) / ei dx-2a^ 2Pko / 

- 2(/ii/?fe + l){2p 2 p k() + a 2 ) / ipie 2 k dx +o(3). 

,/n J 

Thus, the function g(y) in (|4.3p can be written as 

/ ^ P2pk a bi 3 

where &i is as in (|4.10j) . Hence the theorem follows from Theorem A.l in [8]. □ 

n 

When the domain 51 is a rectangle, i.e. O = Yl (0, ij), the b\ in (|4. 10[) for the 

i=i 

Neumann condition can be explicitly expressed in terms of the physical parameters 
fj,i, fj,2, a, and Lj (1 < j ' < n). For example, we consider the case where Q — (0, L). 
The eigenvalues p k and eigenvectors e& of (|3.2I) are given by 

(fc-l) 2 7T 2 (Jfe-lW 

Pfc= ^ , Cfc — cos x, fc=l,2, 



It is clear that fco > 2, and (|4.9|) holds true. We see that 

2 1 ^ 2(feo-l)7T , 1. . 
e ko = 2( 1+COS 1 X > = jp 1 + e 2feo-l]- 

Hence, by (|4.11[) . we have 

(4.16) V = £e x + -qej with j = 2k - 1, 

where 

> - 1 « 2 \ = (--pIpIMpko + 1 

-A -a 2 ) \£ 2 ) \ iplp ko {^Pko + 1 ) 



Ao - (A*iPj + 1) a 2 \ /V\ = (~iplpl {^Pko + 1) 

-A -(P2P 3 +a 2 )J \i] 2 J \ ^p 2 2 p 2 ko {piPk a + I) 

i = — ^— (MiPfco + l)(M2P*o + a2 )- 



PHASE TRANSITIONS FOR THE BRUSSELATOR MODEL 



9 



It is readily to see that 

6 



o, 



" 3 M2Pfc (MlPfeo + 1), 

lAPjPkoiUlPko + 1) 



(4.17) 



??2 = 



a[(M2ft + ct 2 )(piPj + 1) - A*2ftAo] ' 
-fjipko&iPj + l)(MiPfe + 1) 



a[(^ 2 Pj + a 2 ){pipj + 1) - M2ftA ] 
Inserting flTe]) and (|4"T?|) into (|4~10|) , we derive that 



(4.18) 



P2Pk a (PlPk a + l)i 



A i 2Pfc 



-a 



+ 



a 2 {p\Pj + l)p2Pk„ ~ |(MiPfco + !)( 2 M2Pfc + Oi 2 )^2Pj 



Pk 



P.I 



(pipj + l)(p2P 3 + a 2 ) - M2ftA 
x [a 2 (/Xi/) feD + 1) - P2Pk (p2Pk a + a 2 )] 1 , 
(fco - l) 2 ^ 2 



L 2 

4(fc - 1) 2 7T 2 

L 2 



Thus, for the one dimensional domain il — (0,L), the number b\ in (|4. 18|) can be 
equivalently rewritten as 



(4.19) b x = 



P-2Pk + o« - 



P2Pk {^PiP2pj + 4MaPfco - 2a 2 p lPko + a 2 ) 
12pip 2 p 2 o - 3a 2 



x [a 2 (/Jip feo + 1) - p 2 pk (p2Pk„ + a 2 )] 



5. Transition from complex eigenvalues 

As Ai < Ao, the transition of (|2.7[) occurs at A = Ai, and the system bifurcates 
to a periodic solution. 

We first consider the Neumann boundary condition. In this case, Ai = a 2 + 1, 
and we have the following theorem. 



Theorem 5.1. For the problem H2. 3\) with H2.5\) , when Ai < Ao, the transition at 
A = Ai is continuous, and the problem bifurcates on A > Ai to one periodic solution 
which is an attractor. Moreover, the bifurcated periodic solution v x — (v^v^) can 
be expressed as 

v^ = 2[C/3 (X)] 1/2 ^sm(at+^)+o(\f3 \ 1 / 2 ), 
(5.1) 4 

vl = [C/3 (A)] 1/2 a((a - 1) cos at - (a + 1) sin at) + o(\(3 \ 1/2 ), 

where (3 = i(A — and C — (2ira 2 + ^-Ka 4 )^ 1 . 

Proof. We shall verify this theorem by using Theorem A. 3 in [3]. The eigenvalue 
in (|3 .4[) are given by 



Pi W = ^ - A i) ± 2 ^ 4 " 2 (A Al)2 ' ( A i=« 2 + 1 ) 
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Namely, for A near Ai, 



p (\)=Re/3f(\) = ~(\-\ 1 ), 



@f{\i) = ±ia. 
The eigenvectors £ and 77 corresponding to /^(Ai) satisfy 



a 2 






(a 2 - 


h 1) -a 2 ) ( 










a 2 




' m 


(a 2 - 


1- 1) -a 2 ) ( 


V ^ 


£ = 


(^,6) = (« 2 , 


a(l 


7/ = 


(»)i,l2) = (a 2 


, —a 



■'71 

'72 

CI 
6 



It is easy to see that 
(5.2) 

The conjugate eigenvectors £* and 77* satisfy 

' * 2 -(a 2 + l) \ / \ V 



2 -" 2 J \ £2 J " V 7 ?2 

a 2 -(a 2 + l) \ / \ _ ^ / £ 



v 2 



12 / a V £2 



77* = (eta*, 772) = (<*(<* - l),a 2 )- 



leading to 
(5.3) 

It is readily to check that 

(CO = -M*) = 2a 3 f e\dx = 2a 3 \n\, 
(5.4) Jq 

(e,r?*) = (r? ! n=o. 

For the operator G defined by (|2.6I) . we deduce from (|5.2p that for x, y £ K 1 , 
(5.5) 

G(x£ + yrj) 

2Ai 



a 
2Ai 



(x£i +ym) 2 + 2a «i + 2/771) (#£2 + 2/%) + «i + yvi) 2 ( x & + OT2), 

(^6 + Mi) 2 - 2a(a;£i + ym)(x& + ym) ~ «i + OTi) 2 «2 + 27772), 
_ f 2a 3 [(a; + y) 2 + a(a; 2 - t/ 2 )] + a 5 [(a; + y)(x 2 - y 2 ) - a(x + y) 3 }, 
~\- 2a 3 [(x + yf + a(x 2 - y 2 )} - a 5 [(x + y)(x 2 - 17 2 ) - a(x + y) 3 }. 

Because the first eigenvector space Ex = span{£, 77} of (|3.1[) with (|2.5p is invariant 
for the equations (|2.3[) with (|2.5I) . the center manifold function $ vanishes, i.e., 

®(x,y) = 0. 
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Therefore, we derive from (|5.2[) - (|5.5|) that 



(G« + OT + $),£*) 2 , 

; : = ainX + <2l l XV - 


2 ^ 
f- an?y + CL-inX 


+ &2i£ 2 y ■ 


+• ayixy 


(G« + OT + $),//*) 2 
(V,V ) 




+ b 2 ix 2 y 4 


- b\2xy 2 


where 








a 2 o = a(a + 1), 


«02 = 


a(l — a), 




6 20 = a(a + 1), 


002 = 


a(l — a), 




an = 2a, 


hi = 


2a, 






b 03 = 




a), 


a 12 = -Tj^i 1 + 3q 0) 


021 = 


ia 3 (l - 3a). 
2 v ; 



Then, the parameter b in Theorem A. 3 in [3] is 

7T 7T 

o = — (a 20 2 - 020020) + 7- («n«20 + anao2 - 011&20 - &11&02) 
2a 4a 

37T 7T 

+ — (030 + 003) + ^(ai2 + 021) 
3 

= - 7ra 2 (2 + -a 2 ). 

Namely, b < 0. Hence, by Theorem A. 3 in [3], the system (I2.7[) bifurcates from 
(u, A) = (0, Ai) to a periodic solution on A > Ai, which is an attractor. The proof 
of the expression (|5.1[) is classical. Thus the theorem is proved. □ 

Now, we consider the Dirichlet boundary condition. In this case, p\ > and 
Ai = (/ii + ^2) Pi + a 2 + 1. By (|3.5p and (|3.6p it is easy to see that as Ai < Ao we 
have 

M2Pi(M2Pi + a 2 ) < a 2 {piPi + !)• 
Then we define the following parameter 



27ra 2 [/ e?da;] 2 



(5.6) 61 = ^2rr 2L2 fcifr + 1 )^Ip\ + ^Piifiipi + 1) - og) 



Tra 2 f n ejdx 2 
-(2/i 2 pi + 3a ) 



2j n e 2 dx 
2na 2 



(3/xi/Ji + /x 2 /0i + 3)Ai + (/ii/X2Pi + A*2/0i + <*o)-Bi] 



87ra 2 CTn . , o 

"7 K^lPl + M2P1 + 1)^2 - (A*1A*2/>1 + ^2/31 - CT0JS2J 

Jn e i ax 

Ana 2 .. n 2\ x 2/ 

-7 — 2~j [\P1P2P1 + M2P1 ~ ct )A 3 + a {fiipi + p 2 pi + IJ-B3J, 
Jo e i ax 
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where 

cr = [a 2 (^i ) oi + 1) - Hipli^px + a 2 )] 1/2 , 

Ai = (2Ai - (i 2 pi ~ a 2 )B, + a 2 Ci for i = 1, 2, 3, 



fc> 



\ detM kJ Q e 2 dx ' 



^ detM /c det(M^ + 4a 2 ) J n e 2 dx 
R _ y- [Jn e 2 efccM 2 

c = >p (MiPi + l)[/ e 2 e fc dx] 2 

C2 - gdetM fc det(M 2 + 4 CT 2 ) /o e 2 ^^^ + 1 " ~ ^ 



k> -det(M 2 + 4a 2 )J n e 2 dx^^ 
D k = (v2Pk + a 2 ) 2 + 4cr 2 - a 2 {pi + fi 2 )pk - a 2 (a 2 + 1), 
Q k = XiQ 2 - Xi(pi + H2){pk - Pi) - {piPk + 1 - Ai) 2 - 4a 2 . 
Here is the matrix defined by 

-p.ip k + Ai - 1 a 2 

-Ai -M2Pfc - a 2 



(5.7) M fc = 



Theorem 5.2. Let &i 6e i/ie number given by \5. 6]) and Ai < Ao- -for i/ie problem 
\2. 3\) - (2^4\ ), the following assertions hold true: 

(1) The problem undergoes a dynamic transition at A = A 1; which is the Hopf 
bifurcation. 

(2) When b\ < 0, the transition is of Type-I and bifurcates to a stable periodic 
solution on X > Ai, and when b± > the transition is of Type-II and 
bifurcates to an unstable periodic solution on A < Ai . 

(3) The bifurcated periodic solution v x = (^1,^2) can ^ e expressed as 

itf = 2a 2 [- 1 {\)/b l ] 1 l 2 e l ahx(a t + J) + c(| 7 | 1/2 ), 
« 2 A = 2(a 2 + (p 2Pl + a 2 ) 2 )[-2P]i2 ei cos{(Tot + e) + (| 7 |V2 )j 

Ol 

_ a (T + M2P1 + a 2 
= tan -, 

Co - M2P1 - « 
where 7 = (A — Ai)/2. 

Proof. By (|3.3[) the eigenvalues and eigenvectors of p. II) with (|2.4[) at Ai = (/xi + 
A*2)pi + oi 2 + 1 are determined by the matrices M)~ given by (|5.7p . It is clear that 
Mi has a pair of imaginary eigenvalues 

= ±«o" (ct as in (|5.6[l ). 



(5.8) 
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Let £, f] e M 2 be the eigenvectors of M\ satisfying 

Mi| = <T fj, Mifj = -<ro£. 

Then, by (|3 .3[) the eigenvectors of (13. If) corresponding to /^(Ai) are given by £ 
£ei and 77 — fje±. It is readily to check that 

(5.9) £ = (£1,62) = (a 2 ei, (cr - M2P1 - a 2 )ei), 

(5.10) 77 = (771,772) = (a 2 ei, -((T + /j 2 Pi + a 2 )ei). 

We consider the conjugate eigenvectors £* = £*ei and 77* = 77* e\ with 
M?t=<roV*, M*fj* = -a t, 
where M* is the transpose of Mi . Direct calculation shows that 

(5.11) C = = ((<7 + M2P1 + « 2 )ei,a 2 ei), 

(5.12) 77* = (771,772) = ((-co + /U2P1 + a 2 )ei,a 2 ei). 



It is easy to see that 
(5.13) 



(CD = -fof?*) = 2a 2 cr / n e 2 da;. 



Let u = x£ + yr/ + $(x, y) € i? be a solution of (f2~3]) - (|2~4l) at A = Ai, and $ be 
the center manifold function. By (|5.13[) the reduced equations of (|2.3p - ()2.4j) read 

dx 1 

(5.14) * > 
at (77,77 ) 

where the operator G is given by 

(5.15) G(u) = G 2 (u)+G 3 (u), 
and Gk {k = 2, 3) is a fc-multilinear operator defined by 



G 2 {u, v) = 2 ( —uivi + auiv 2 ,- ( — iti«i + cwitj 2 

(5.16) G 3 (m, tj, io) = (uitjiw 2 , -UiViW 2 ), 

G 2 (u) = G 2 (u,u), 
G 3 (u) = G 3 (u,u,u). 



)) 
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Based on ([5T9 ]) - (|533]l and (f57T5^) - ff5TT6]) . (|5~l"4l) are rewritten as 



(5.17) 



where 



dx 
~dt 



dt 



= - &oy + fl20^ 2 + anxy + a 02 y 2 + a 30 x 3 + a 21 x 2 y 
+ a 12 xy 2 + a 03 y 3 + j^- (G 2 (C <f ) + G 2 ($, 0, O 
+ ^py (G 3 (t?, <*>) + G($, r?), ,,*) + o(3), 

=CT a; + & 20 x 2 + 6n^2/ + b 2 2V 2 + hox 3 + b 2 ix 2 y 

+ b 12 xy 2 + b Q3 y 3 + 7 _£--(G 2 (£, $) + G 2 ($, CW) 
(7?, 77*) 

-(G 2 (77,$) + G 2 ($,77),^)+o(3), 



(r/,7?*) 



(G 2 (CDO 4^ 

(CD a aJ n eidx 

«n = = ? — 3-7- ^0 + M2P1 (M1P1 + 1 , 

(CD cr J o efrfa; 

(G 2 (T?,T?),D , , w x 

002 = = j — rr^ + M2P1 (M1P1 + 1 - cr , 

(G 3 (CCDO °? In 4dx , ^ 2 
fl 3o = 77— ttn = o r 2^ g o + M2P1 (o-o - M2P1 - or), 

«2i =7^ry (G 3 (C C r?) + G 3 (C 77, + (G 3 (rj, C 0, D 
-(<t + M2Pi)(^ - 3^ 2 Pi - 3a ), 



2c7o Jn e f^ 
2cto e\dx 



«i2 = 7F ^-r(G 3 (C »?, r?) + G 3 ( V , C r?) + G 3 ( v , V , £), D 

(CD 

(cr + M2Pi)(ct + 3^ 2 Pi + 3a 2 ), 



a °3 = tft^; = _ o — r ia (°"° + M2Pi)(°"o + M2P1 + or), 

(CD 2lJ oJ n efdx 
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and 



b 20 = ; rr = r 2 , (go - M2P1XM1P1 + 1 + °"o), 

(77,77*) ao J n e(dx 

(G 2 (ti,ti),t)*) oj^efrfx 
&02 = (77,77*) = ao/nefdz^ " + 1 " 

(77,77*) 2a J n efdx 

621 =^7(03^, £, r?) + G 3 (£, 77, + (7(77, & 0: »?*) 
(77,77*) 

a 2 f efdx 2 
-(o-o - M2Pi)(o-o - 3/x 2 pi - 3a ), 



2ct / n e 2 cte 

' -(G 3 (t7, 77, + 63(77, e, r?) + G 3 (£, 77, 77), 77*) 



■' ! ' 1 -(ct - p 2 pi)(cro + 3/^ 2 pi + 3a 2 ), 



^03 



2cr o / n e 2 dx 
_(G 3 (77, 77,?7),77*) a 2 Lefdx 



(cr - fX2pi)((To + P2P1 + a 2 



(77,77*) 2ct J n efdx 

We are now in a position to derive the center manifold function $. By (A. 10) in 

(5.18) $ = $ 1 + $ 2 + $3 + o(2), 

where 

- L Al $i = P 2 [G 2 (e, 0^ 2 + (G 2 (f, r?) + G 2 (7?, + G 2 (t7, 77)7/] , 

- [L 2 Al + 4^]L Al $2 = 2a 2 P 2 [(G 2 (£, - G 2 (t7, 77))( y 2 - a; 2 ) - 2(G 2 (£, 77) + 62(77, £))xy] . 
[L\ x + 4a 2 ]<J> 3 = a P 2 [(G 2 (£, v) + G 3 (tj, OX?/ 2 - ^) + 2(G 2 (£, - G 2 (t7, r,))xy] , 



L\ 1 = A + B\ 1 is the linear operator defined by (|2.6p . P 2 '■ H ^ E 2 the canonical 
projection, and £2 = {" € = 0, (u, 77*) = 0} is the complement of £1 = 

span{£, 77} in H. Note that the eigenvectors of L\ x satisfy 



<pk = 4>kei, 4>k € M 2 , fc = l,2,---, 

= /?/c0fc (-Wfc the matrix as in (|5.7|l ). 
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Hence, we obtain from (|5T9 ]> . (jSTTOjl . f[57T6)) and ([538]) that 

$1 =2a 3 [(p, 1 p 1 + a + l)x 2 + 2(/iipi + X)xy + (pipi - a + l)y 2 



k>i /n 



0> 2 =8a 3 a 2 la (y 2 - x 2 ) - 2(/x lPl + V J " 1 {-M k )-\M 2 k + 4a 2 ) , 

Jn e fe a;E V 



efe, 



$ 3 =4a 3 a [(p 1 p 1 + l)(y 2 - x 2 ) + 2<j xy] V Jsl * (M 2 + 4a 2 )- 1 _ | < 

fc> i Jfi 6 *^ 

Direct calculation shows that 
1 



(-M fe ) 



-1 



detM, 



Ik 



M2Pfe + a 2 a 2 



-Ai ^iPfc + 1 - A 



v k 0J det(M 2 + 4(7 2 ) 
Thus we have 



(p 2 Pk + a 2 ) 2 + 4cr 2 - Aia 2 a 2 (pi + p 2 )(Pfc - Pi) 

-A1O1 + M 2 )(p/c - Pi) (piPfc + 1 - A1) 2 + 4er 2 - Aia 2 



$1 = 2a 3 [(pipi + 1 + (T )x 2 + 2(pipi + l)xy + (nipi + 1 - cr )y 2 
$ 2 = 8a 3 cr 2 [-cr a; 2 - 2{ji l p l + l)xy + <r y 2 } ( ^ 
$ 3 = 4orV [-(AiiPi + 1)-^ 2 + 2(T -Ty + (piPi + l)y 2 } 



£1 



E 3 
F 3 



where 



rp _ In e i e kdx 

El -h^M k f n e 2 dx^ 



OC 



^ = -gdetM fc/ ,e^^ + 1 ^ 

p - Jo ejekdx 

h3 ^et{Ml + Aai)S n eldx Dkek > 

^-^ i det(M 2 +4a 2 )/ n e 2 dx^ efc - 
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Inserting $ = $1 + <f> 2 + $3 + o(2) into (|5.17[) we derive that 

~[£ = ~ (Toy+ X] a ij x l y j + ^2 a kr x k y r + o(3), 

' 2<i+3<3 fe+r=3 

( 5 - 19 ) dy 

-^ = cr a;+ b i jx l y j + ^ b kr x k y r + o(3), 

2<i+j<3 fc+r=3 

where a.y and 6^ (0 <i,j< 3) are as in (|5.17|) . and 

030 = (o~o + ^Pi) + 1-1- o-q)^ - 40^/2 - 2ct (mi/°i + 1)^3], 

do J n ejdx 

ai2 = 2a r + 2 ^ 2Pl) [(/ii/Qi + 1 - <7 )Ii + 4<T^2 + 2<T (/x 1 pi + l)/ 3 

+2(mi/0i + 1) Ji - Sai(m Pl + 1) J 2 + 4(t 2 J 3 ], 

603 = 2Q ~ 2 M , 2Pl) [(^iPi + 1 - <ro)Ji + + 2a Q (^ lPl + 1)J 3 ], 

ct J n ejdx 

b 21 = 2^Z^)[ 2 ( WP1 + 1)/! - 8a 2 ( MlPl + l)/ 2 +4(T 2 / 3 
cto J e^x 

+(/Xi/>i + 1 + <7 )Ji - Acrl J 2 - 2cr (MiPi + 1) J 3 ], 

where 



Jj = A; + (T - / Eie\dx, 
Jn 

Ji = Ai — cr / Eiejdx, 
Jn 

Ai = (2Ai - /i 2/ 9i - a 2 ) / £ j; e 2 cfa; + a 2 / F { e\dx, « = 1,2,3. 
Jo Jsi 

Then, the number 

7T 7T 

&i =- — (002^02 - 020^20) + 1 — (ana 2 o + «ii«02 - &11&20 - ^11^02) 
2ao 4cr 

+ ^-(a 3 o + ^30 + 030 + ^30) + ^( fl 12 + &21 + ai2 + &21) 

is the same as in (|5.6[) . Thus Assertions (l)-(2) of this theorem follows from Theo- 
rem A. 6 in [TP] , 

It is known that the bifurcated periodic solution near A = Ai takes the form 
(5.20) v x = x(t)£ + y(t)r, + o(\x\ + \y\), 

where £, 77 are as in (|5.9p and (|5.10() . and (x(t),y(t)) is the solution of the following 
equation 

1 

— = j(X)x - a (X)y + j—-(G(x^ + yi]\ + $ A ), ft), 

-| = (t (A).t + 7 (A)y + 1 + OTa + $a), r? A ), 

where £,\,r]\ are eigenvectors of L\ corresponding to the first complex eigenvalues 
fif = 7±zcr , and £ A ,r? A ^ ne conjugate eigenvectors. This solution (x(t),y(t)) near 
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Ai is of the form 
(5.21) 



x{t) = 

y(t) = 



7(A) 

7(A) 
bi 



1/2 



1/2 



cosfToi + o^MI 1 / 2 ), 
sina i + O (| 7 (A)| 1 /2), 



where b\ is as in (|5.6j) . Therefore, Assertion (3) follows from ()5.20|) and (|5.21| . The 
proof is complete. □ 



6. One-dimensional case 

When the containers f2 are taken as rectangles, the criteria in Theorems 14.21 and 
I can be simplified. For simplicity, we consider here only the one-dimensional 
case: 

(6.1) SI = (0, L) C R 1 . 

The eigenvalues pk and corresponding eigenvectors of 
( k 2 n 2 /L 2 for b.c. (HU) , 



(6.2) 
(6.3) 



pk = 



(fc - l) 2 7r 2 /i 2 for b.c. (1231) . 
. k?tx 

(k — l)nx 



are given by 
(fc = 1,2, - ■ -), 



sin ■ 



cos • 



for b.c. ([231) 
for b.c. (|2~B1 



Thus the two critical numbers Aq and Ai in 



and (|3.6p are given by 



(6.4) 
(6.5) 



A = min fc 2 



Ai 



\jL\k 2 -K 2 
L 2 

1 1- M2) 
1 



a 2 L 2 



Ml" 



1 



p 2 k 2 TT 2 [l 2 

1 ' 1 for b.c. (EH), 
for b.c. (I2~5l). 



It is known that the criterion b\ in Theorem 14.21 is valid only for the free boundary 
condition, which can be expressed explicitly by (|4.19l) . Likewise, for the number 
defined by (|5.6I) we have the following explicit expression 



(6.6) 
with 

(6.7) b 



b\ = -^j—°0, 



4 3 x L 



9<J 2 1T 2 



{piPl + l)(M 2 Pl + IpiPlPl + 2 M2Pl ~ CT ) 



3i 
16 



(2^ 2 Pi + 3a 2 ) + 2{3fi lPl + M2P1 + 3)Ai 



-Sa^ipipi + p 2 pi + l)A 2 - A(p,ip, 2 pl + P2P1 - o- 2 )A 3 
+2(pip 2 pl + V2P1 + u a )Bi + 8o- 2 (/ii/i 2 /9 2 + p 2 pi - &o)B 2 
-4(7o(/iipi + p 2 pi + 1)5 3 , 
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where A, = (2/xipi + p 2 pi + a 2 + 2)B t + a 2 C t (i = 1, 2, 3), and 



B 1 = 



/'2 



k=l 



detM 2fc+1 [(2fc + l) 2 -4] 2 L 2 ' 



R =W"V (fJ-2P2k+i + a 2 )D2k+i + a 2 Q2k+i 

2 fci detM 2fc+1 det(M 2 2 fe+1 + 4 ( 7 2 ) 7 r 2 (2fc + l) 2 [(2fc + l) 2 - 4] 2 ' 



— det(M 2 k+1 + Aa 2 )Tr 2 (2k + l) 2 [(2k + l) 2 - 4] 2 : 
MiPi + 1 



C 2 = -32L^] 



^ detAf 2fc+1 7r 2 (2fc + l) 2 [(2fc + l) 2 - 4] 2 ' 

(piPi + P2P1 + a 2 + 1, -D 2 fc + i - {piPk - pipi - P2P1 - a 2 )Q 2 k+i 



detM 2k+1 det(M 2 k+1 + Aa 2 )n 2 (2k + l) 2 [(2fc + l) 2 - 4] 2 ] 



C 3 = 32LJ2 



k=l 

n 

W2k+l 



k ^ det(M 2 k+1 + 4a 2 )n 2 (2k + l) 2 [(2fc + l) 2 - 4] 2 ' 



(. b -^ A o - — + .. ,„2^ 2 + — a + L 



Let Ao in (|6.4p achieves its minimum at the integer /cq, i.e. 

p,\k 2 -K 2 a 2 L 2 pi 2 
L 2 p 2 k 2 Tr 2 p 2 
Then fco > 1 satisfies 

ah 2 

(6.9) fc (fc - 1) < , - < fc (fc + 1). 

T^\JP\P2 

To see this, note that the function 

., . p\-K 2 x a 2 L 2 ui 2 
X(x) = r9 + — + —a 2 + 1 

L 2 - p 2 TT Z X p 2 

has its minimum at xq = aL 2 /ir 2 y?pip 2 , and 

d\ J < if x < i , 
da; \ > if x > xq. 

It follows that either kg = to or ko = m + 1, such that m = y/xo ~ e for some 
< e < 1; namely 

m 2 < x < (to + l) 2 . 

It follows that 

{to if A(m 2 ) < A((m + l) 2 ), 

to and to -)- 1 if A(m 2 ) = A((m + l) 2 ), 
to + 1 if A(m 2 ) > A((m+ l) 2 ). 

We infer from (16. 10)) that 

ah 2 



aL 2 



< m(m + 1) => ko = m, 
> to(to + 1) =>■ ko = m + 1, 
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which yield the inequalities (|6.9I) . 

In the following, we compare Ao with Ai in terms of the parameters fn, fJ>2,cx, 
and L. We proceed in two cases. 



The case where > H2- Then, from (|3.5|) - (|3.6p , we see that 
(6.11) Ai < A V£, a > for b.c. flUg) . 

For b.c. (|2.4I) . we can prove that 

A > Ai if < L < L c , 

Ai < A , if L > L c , 

where 

(6.13) Ll = ^-[yj (hi - H2?c? + ij4 - afa - (Mi)]. 

In fact, from Aq = Ai, we derive the critical scale L c as 



Ptt 2 
r2 _ li o 7r 



(fii - H2) 2 a 2 + -rr (Mi + M2 - feoMi) ~ a (Mi - Ma) 



^o 



provided that 



Mi + M2 > fcoMi- 



Since fii > M2, it implies that ko = 1. Hence L c is as in (|6.13[) . and ()6. 12|) holds 
true. 

We remark that by (|6.9p . L 2 in (16.131) has to satisfy the inequality 

L 2 < 2ir 2 y/HiH2/u, 

which holds true for hi > M2- 

The case where hi < ^2- For b.c. (|2.4p . we deduce from A = Ai the following 
two critical scales: 



k 2 ir 2 



(M2 - Mi)« T W(M2 - Mi) 2 " 2 



4/x 2 

1,2 



O0M1 - Mi - M2) 



with L 2 x < L 2 2 

It is easy to see that for the boundary condition (|2.4j) . 

•15) 

' k 2 a 2 (H2 ~ Mi) 2 < 4/i 2 (fc 2 Mi - Mi - M2), 
or 



Ai < Ao if < 



(6.16) 
A < Ai 



if 



fc 2 a 2 (/i 2 - Hi) 2 > 4^2(fc 2 Mi _ Ml - M2) and 
< L 2 < L 2 C1 or L 2 C2 < i 2 , 



L 2 C1 < L 2 < L 2 C2 for < 4^ 2 (fc 2 ^i - H1V2) < k 2 a 2 (p2 - Mi), 
or < L 2 < L 2 for k 2 Hi < Hi + M2- 



For the boundary condition (I2.5[) . we obtain two critical scales as 



(6.17) 



2a 



(M2 - Mi)" T V (M2 ~ Mi) 2 " 2 - 47*1 



M2 
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such that 

!(M2 - Ml) 2 " 2 < 4^1/^2, 
or < L 2 < l\, or l\ < L 2 , 
for (p 2 ~ Mi) 2 " 2 > V1M2, 
(6.19) A < Ax for b.c. flU} if < L 2 < l\ for (/x 2 - Mi) 2 a 2 > 4/iiju 2 . 

7. Physical Remarks 

We now discuss the phase transition of Brusselator by using Theorem 14. 11 5 . 2 1 for 
the one-dimensional case (|6.ip . 

Dirichlet Boundary Condition. When /^i > fi 2 , by (|6.12p . the system 
(|2.3[) - (12.4[) has a transition to steady states provided < L < L c , and to periodic 
solutions provided L c < L. 

Physical Conclusion 7.1. Let \i\ > /i2- Then for the system 112. 3\) - [2~4\ ), we have 
the following conclusions: 

(1) When < L < L c , the transition of 12. Sty -{2.$ at A = A is of Type-Ill, 
and there is a saddle-node bifurcation at some < A* < Ao- In other wards, 
the basic state u$ = (a, A/a) is stable for < A < A* 7 is metastable for 
A* < A < Ao, and is unstable for A* < A. Moreover, if X* < X, there are at 
least two metastable equilibrium states. 

(2) When L c < L, this system undergoes a dynamic transition at X = Ai 
to periodic solutions. In particular, there exists an Lq > L c such that if 
L c < L < Lq, the transition is of Type-II, and there is a singular separation 
of periodic solutions at some X < X\. If Lq < L, the transition is of Type-I, 
and the system bifurcates from Uo to a stable periodic solution on X > Ai, 
which is expressed as u\ = (u^'u^) with 

u$ = a + 2a 2 ^^-p- sin ^ sin(cr i + J) + o(|A - A^ 1 / 2 ), 

u\ = ~ + 2( CTo 2 + (a 2 + ^) 2 )^^isin ^ cos(a i + 0) + o(\X X^ 2 ). 

This periodic solution provides a spatial-temporal oscillation of the Brusse- 
lator. 

The first conclusion is due to Theorem 14. 1\ the existence of global attractors, 
and the fact that Uo = (0, 0) is a unique steady state solution of (|2.3[) - (|2.4[) at A = 0. 

The second conclusion is based on Theorem 15.21 and the following analysis on 
the criterion bo given by (|6.7[) . We know that 

Ai — A — > as L — > L c , 

which implies 

(To — > as L — > L c . 
It follows from (16.71) that 60 +°° f° r L L c + 0. Therefore 
(7.1) b a > VL C <L< L , 
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for some Lo > L c . On the other hand, pu — > (L — > oo). Hence, when L — > oo, 

(7.2) b -»■ - ^ - + 2(3a 2 + 6 + c7o)Bi - 8a 2 (a 2 + a + 2)B 2 

+ Aa 2 {a 2 + 1)B 3 + 6a 2 C\ - 8a 2 a 2 C 2 + 4a 4 C 3 . 

Note that 

cr — > a, D k — > 3a 2 , Q fc -> —3a , 

detM fe -> a 2 , det(M^ + 4a 2 ) -> 9a 4 , Bi -> 0, 

3a z a z 

for L — ¥ oo, where 

oo 

g _ i 

^7T 2 (2fc + l) 2 [(2fc + l) 2 -4] 2, 

Thus, in view of lfT2]) . 

64 9a 2 
l 9tt2 + ~L6 

Hence 



6 -» -£[7^2 + "77" + 64£; ] < as L ^ oo. 



(7.3) b < VLi < L, 
for some L\ > L . 

From the physical point of view, it is reasonable to consider the case where b 
changes its sign only once in (L c ,oo). Hence, physically, we have Lo = L±. Thus, 
we derive from (|7.1[) and (|7.3p the second conclusion. 

Now, we consider the case where /ii < fi 2 by the following two examples. We 
take 

(7.4) hi = 2 X 1(T 3 , /i 2 =4x 1(T 3 . 

Example 7.1. Let (f7T4"|) hold true, and a = 2, L — 4. Then we obtain from (RUty 
that fco — 34, and 

k 2 a 2 (fi 2 - Mi) 2 < 4/x 2 (fco/Xi - pi - A*2)- 

In view of (|6. 15|) . 

Ai < A for b.c. (gH) . 
The number 6o in (|6.7p is given by 
/ 64 9\ 

b ^ -( — + -) L + 40Bi - 256B 2 + 80B 3 + 24Ci - 128C* 2 + 64C 3 , 



and 



9tt 2 



2 11 

Si = — x 10~ 3 L, B 2 ^ x 1CT 3 L, 

25 36 x 25 

J 27 x 25tt 2 ' 9 x 25tt 2 ' 

Co — ~t L . Cr> = — L . 

2 27 x 25tt 2 ' 3 27 x 25tt 2 
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It is easy to see that 60 < 0. Then, by Theorem 15.21 the phase transition of 
(|2.3p ~ (|2.4[) is of Type-I, and this system undergoes a spatial-temporal oscillation on 
A>Ai. 



Example 7.2. Assume (|7.4|) and a = 3, L = 4. Then k = 41, and 

k%a 2 (n 2 - Mi) 2 > 4/^2 (&oMi - Mi - M2) > 0. 
Thus, we can obtain two critical scales L\ and L 2 2 in (|6. 14|) as follows: 



Hence 



L 2 C1 = 11.07, L 2 C2 = 22.14. 



L\ < L 2 = 16 < L 2 C2 , 



which leads (by (|6.16|) ') to 

(7.5) A < Ai for b.c. (El 

It is clear that 



~l o dx = I sin 3 -^—dx ^ 0. 



si Jo 



L 



By ([73]) and Theorem H7TJ the system ([23)1 and (|2.4p has a Type-Ill transition at 
A = Ao, and there is a saddle-node bifurcation at some A*(0 < A* < Ao). 

Neumann Boundary Condition. Consider the case where /zi > fi 2 - Based 
on (|6.15p and Theorem 15. 1[ we have the following physical conclusion. 



Physical Conclusion 7.2. Let [i\ > M2- Then the system \2. 3\) with \2. 5\) has a 
Type-I transition to periodic solutions at A = Ai, i.e., a spatial-temporal oscillation 
occurs in the Brusselator for A > Ai . 

For the case where /ii < M2, we have the following example. 

Example 7.3. Under the same conditions as in Example 1 7. 2 [ ko — 41 and the two 
critical scales l\ and l 2 2 in (|6.17p are given by 

l 2 Cl = 11.06, 1% = 22.12. 

Hence, l\ < I 2 = 16 < l 2 2 , which implies, by (jnHSJ), that 

(7.6) A < Ai for b.c. ((2~5l) . 

On the other hand, by (|6.2[) and (|6.3p . we have 

L2 

_ v o 

Pk + 1 



|^-10 3 , J el 0+x dx = J cos 3 ^j-dx = 0. 



in Jo 

Thus, it is easy to check that the number b\ in (|4.10| . which is also given by (|4.19[) . 
is negative, i.e., 

(7.7) 61 < in (j4T0)l . 

By (f7T6]) - (T777|) and Theorem [431 the system (f2T3]> with (|2~5|) bifurcates on A > A 
to two stable steady states v± as given by ()4.12j) . It shows that the Brusselator 
undergoes a transition at Aq = 9.8. 
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